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Abstract:- In this paper we study the existence of a periodic solution for nonlinear system of differential 
equations with pulse action of parameters. The numerical-analytic method has been used to study the periodic 
solutions of the nonlinear ordinary differential equations that were introduced by Somioleko And the result of 
this study which is the using the periodic solutions on a wide range in difference processes in industry and 

technology. 



I. INTRODUCTION 

There are many subjects in physics and technology using mathematical methods that depends on the 
nonlinear differential equations, and it became clear that the existence of the periodic solutions and it's 
algorithm structure from more important problems in the present time. Where many of studies and researches 
dedicates for treatment the autonomous and non-autonomous periodic systems and specially with the integral 
equations and differential equations and the linear and nonlinear differential and which is dealing in general 
shape with the problems about periodic solutions theory and the modern methods in its quality treatment for the 
periodic differential equations. 

Somioleko [6] assumes the numerical analytic method to study the periodic solutions for the ordinary 
differential equations and its algorithm structure and this method include uniformly sequences of the periodic 
functions and the result of that study is the using of the periodic solutions on a wide range for example see [4, 5, 
6]. 

Consider the following system of nonlinear differential equation, which has the form: 



dx 

— = Ax + f(t,x,y) , t±ti , Ax 
at 



dy 

— = fix + g(t,x,y) , t±ti , Ay 



(1) 



Where x £ D x £ R n , y £ Dp £ R n , D^is a closed and bounded domain. 

The vector functions f(t,x,y) and g(t,x,y)are defined on the domain: 

(t,x,y) £ R 1 x D x x D p = (-00,00) x D A x Dp . . . (2) 

Which are continuous int , x , yand periodic in t of period T, where Dp is bounded domains subset of 
Euclidean spaces R m , and the functions l t {x,y) , Gj(x,y)are continuous in the domain (2), where li+ v {x,y) — 
h( x >y)> G i+p (x,y) = Gj(x,y) and t i+p (x,y) — t t +7 for p is a positive integer and {t;} is finite positive 
sequence of numbers. 

Suppose that the vector functions in (l)are satisfying the following inequalities: 

, max \\f(t,x,y)\\<M 1 , max \\g(t,x,y)\\ < M 2 . . . (3) 

te[0j] te[0j] 

max ||/;(x,y)|| < M 3 , max \\G i {x,y)\\<M A ... (4) 

(xjOeD^x Dp (x,y)€D A x Dp 

\\f{tX,yi)-fit,x 2 ,y 2 )\\<K 1 \\x 1 -x 2 t+K 2 \\y 1 -y 2 \\ ... (5) 
Il5(t,x 1 ,y 1 )-,g(t,x 2 ,y 2 )|| < - x 2 \\ + L 2 \\ yi - y 2 \\ ... (6) 

ll/j(*i.yi)-U*2.y 2 )ll <ff3ll*i-*2ll + ff4llyi-y2ll ■ ■ ■ (7) 

l|G«(*i.yi)-G<(*2.yi)ll <M*i ~x 2 \\ +L 4 \\ yi -y 2 || ...(8) 
Where t £ R' , x ,x 1 ,x 2 £ D x , y ,y 1 ,y 2 £ Z^and M t ,M 2 ,M 3 ,M 4 ,K 1 ,K 2 ,K 3 ,K 4 , 
Li ,L 2 ,L 3 ,L 4 are apositive constant ,|| . ||= max^^ 1. 1 . 

Let, f> are a positive parameter defined in (2), continuous and periodic at t , s , t and satisfy both following 
inequalities: 
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||e* (t-s) || <H 

|| e /?(t-5)|| < F 

Where H and F are a positive constants. 
We define the non-empty sets as follows: 



... (9) 



// XTH 2 T\ ( XT \\ \ 

(1TF 2 T\ ( (3T \\ 

-rr^j M2+pF ( 1 + iT^)j . 



% = D p - | | FT 



. . . (10) 



Furthermore, we suppose that the larges Eigen-value for the following matrix: 
Ao = 

W 3 W 4 . 
is less than one, i.e. : 



Qmax 

where 



W t + W 4 + 4Wi + W4) 2 + 4(W 2 W 3 - W1W4) 



■ ■ ■ (11) 



w 2 



I XTH 2 T\ i XT \ 

l XTH 2 T\ 1 XT \ 

= ( Hr "iT^j^ + 2p// ( 1 + iT^)^ ' 

/ (3TF 2 T\ 1 (3T \ 

w ^{ FT -iT^) L ^ 2pF { 1 + YT^)^ • 

( RTF 2 T\ 1 RT \ 

/ XTH 2 T\ ( XT \ 

M s = {"T- T ^ W )M 1 + 2 P h(i + TT7W )m 3 , 

( (3TF 2 T\ 1 PT \ 

M6 = r-r^r 2+2pF ( 1+ ^) 



Approximation solution of (1) 

The investigation of approximation solution of the system (1) will be introduced by the following 

theorem: 
Theorem 1 

If the system of nonlinear differential equations with pulse action (1) satisfy the inequalities (3) — (8) and the 

conditions (9) , (10) then the sequences of functions : 

t 

* m +i(t.*o.yo) = *o + j e 1(t_s) [/(s,x m (s,x 0 ,yj,y m (s,x 0 ,y 0 )) - 
0 

T 

r ^ 

~ I l + e XT eA(T ~ s) f( s ' x rn( s ' x o,yo),y m (s,x 0 ,y 0 ))ds]ds + 

0 

+ ^ e 1(r - s) / i (x m (t,x 0 ,y 0 ),y m (t,x 0 ,y 0 )) - 

0<t,<t 

~ 1 + e AT ^_ j eUT ~ S)l i( X m( t ' X o>yo)>ym(t,X 0 ,y o )) . . . (12) 



i=l 

with 

x 0 (t,x 0 ) — x 0 e At , m =0 , 1 , 2 , 
end 

t 



y m +i(t,*o,y 0 ) = Vo + j e p[t s) Lf(5,^ m (s,x 0 ,y 0 ),y m (5,x 0 ,y 0 )) - 
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T 

p 



r p 

- J 1 ~ efST e p(T ~ s) f{s,x m {s,x 0 ,y 0 ),y m {s,x 0 ,y 0 ))ds]ds + 
o 

+ ^ e^ T ~^G i (x m (t,x 0 ,y 0 ),y m (t,x 0 ,y 0 )) - 

0<tj<t 

p 

- t ~ e p T ^ ^"^ G i ( X m (*> X o> Vo). Vm (f, X 0 , y Q )) . . . (13) 



With 

y 0 (t,x 0 ) = y 0 e pt , m=0,l ,2, , 

Are periodic in t of period T, and are uniformly convergent as m -» oo in the domain : 

(t,x 0 ,y 0 ) £ R'x% xD pf ... (14) 

To the limit function ^(t, x 0 ,y 0 )andy a3 (t, x 0 ,y 0 ) define in the domain (14), which is periodic in t of period T 

and satisfying the system of integral equations : 
t 

x(t,x 0 ,y 0 )= x 0 + j e Ut - s) [f(s,x(s,x 0 ,y 0 ),y(s,x 0 ,y 0 )) - 
o 

T 

r ^ 

~ I l + e 1T eA(T ~ s) f( s ' x ( s ' x o>yo),y(s,x 0 ,y 0 ))ds]ds + 
o 

+ ^ e A(r - s) /;(x(t,x 0 ,y 0 ),y(t,x 0 ,y 0 ))- 

0<tj<t 

i=l 

and 

t 

y(t.^o.y 0 )= y 0 + Je^ (t - s) [/(s,x(s,x 0 ,yJ,y(s,x 0 ,y 0 )) - 
o 

r 

~ I l + e /?r e ^ T ~ s) f(. s ' x ( s ' x o,yo),y(s,x 0 ,y 0 ))ds]ds + 
o 

+ ^ e^ T - s) Gi{x{t,x 0 ,y 0 ),y{t,x 0 ,y 0 ))- 

0<t ( <t 

- Yf^2^ enT ~ s)G ^ x{t ' Xo ' yo) ' y{t ' Xo ' y °^ 1 ■ ■ ■ (16) 

i=l 

Which is are unique solutions of the system (1), provided that : 

/ ATH 2 T\ / AT \ 

ll*°(t,x 0 ,y 0 )-x 0 || <(//T- T - r ^ F jM 1 + 2p//^ + T -^ F jM 3 . . . (17) 

/ RTF 2 T\ I RT \ 

\\y°(t,x 0 ,y 0 ) - y 0 || < ( FT - M 2 + 2pF (l + y^w) M 4 ■ ■ ■ ( 18 ) 



< A^f-A,,)" 1 ^ , ... (19) 



,\\ya,(t,x 0 ,y 0 ) - y m (t,x 0 ,y 0 ) 
for all t £ [0, T] ,x 0 E D xf , y 0 £ D pf , when: 
^(t) W 3 (t) N 



A 0 = 



y JV 2 (t) iV 4 (t) 
And where 

ATH 2 T 

N ' 



( ATH T\ i AT \ 

^) = (hT- t ^ w )k 1+ 2 P h(i + i ^ w )k 3 , H^max^G) 

/ ATH 2 T\ ( AT \ 

N 2V = ["T-Y^) K2 + 2pH \ 1 + T^) K * ' W2 ^ t f^ m) 
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/ /3TF 2 T\ ( /3T \ 

/ RTF 2 T\ I RT \ 

// ATH 2 T\ i AT \ \ 

{ HT -i^) M ^ 2pH ( 1+ r^) M A 



Wo = max JV 3 (t) 
s te[0,r] dv 

W 4 = max JV 4 (t) 
4 te[0,r] 4V 



F7 



RTF 2 T\ / RT \ 

+ 2pF ( 1 + IT^) M V 



1 + 



M 7 



Proof: 



Setting m = 0 and using (12) , and the condition (9), we get 
t 

||x 1 (t f * 0f y 0 )-* 0 || = |||e 1(t - s) ||||/(5,x 0 ,y 0 )||d5- 
o 

T 

- r ^ If |||e 21(r - s) ||||/(5,x 0 ,y 0 )||d5+ ^ || e 1(r - s) ||||/ ; (x 0 ,y 0 )||- 

0 0<tj<t 

At P 

i=l 

II At// 2 || r At// 2 r 

< - 1 + e xr I H/C 5 ^o<yo)ll^- 1 + e/lr I ll/(s,x 0 ,y 0 )||ds + 

0 0 

ZAtH \ 1 
IIA(*o.yo)ll -i^irZ. l|/ ' (Xo,3 ' o)l1 



0<tj<t 



Hence 

/ ATH 2 T\ ( AT \ 

||% 1 (t,x 0 ,y 0 )-x 0 || <(//T- r ^jM 1 +2p//^ + T - r ^rjM 3 ... (19) 

So thatXj (t, x D , y 0 ) 6 D x , for all t£fl',i 0 £ D^.and by mathematic induction we get: 

/ ATH 2 T\ I AT \ 

\\x m (t,x 0 ,y 0 ) -x 0 \\< \HT - Y^F) M i + 2 ? H + Y+^J Ms 

and 

/ RTF 2 T\ { RT \ 

\\yi(t,x 0 ,y 0 )-y 0 \\ <\FT-^^\M 2 +2 V F[l+^-^jM, . . . (20) 

Hence (t, x 0 , y 0 ) 6 Dp , for all t £ R , y 0 6 Dpf . and by mathematic induction we get : 
/ RTF 2 T\ ( RT \ 

lly m (t.* 0 ,y 0 ) - ydl < (ft - m 2 + 2 V f (i + z^f) m 4 

Then x m (t,x 0 ,y 0 ) £ D A ,y m {t,x 0 ,y 0 ) £ Dp , x 0 £ D Xf ,y D £ D pf . 

We claim that the sequence of functions (12) and (13) are uniformly convergent on the domain (14). 

By using (19), and when m = 1 , we get 

t 

ll*2(t.*o.y 0 )-*i(t.*o.yo)ll = j"||e 1(t - s) ||ll/(X*i,yi)ll^- 

0 

t T 
0 0 

At P 

+ £ ||e 1(r - s) ||H/ ; ^i,yi)ll+ IT ^rZll eA(r " S) ll ||/i(Xl ' yi)l1 " 

0<tj<t 1=1 

t t 
- j ||e 1(t - s) ||||/(5,x 0 ,y 0 )||d 5+I -^ If j ||e 1(r - s) ||||/(5,x 0 ,y 0 )||ds + 



5 + 
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+ Y^fjWe^WWfUs.x^Wds- £ ||e 1(r - s) ||||/ ; (x 0 ,y 0 )||- 

0 0<tj<t 

p 

- IT^Zll el(r " s) llll«^yc)ll , 

i=i 

/ ATH 2 T\ i AT \ 1 

{ HT -TT^) K i + 2pH { 1 + TT^) K *\ 

ATH 2 T\ i AT \ 1 

^-^)K 2+ 2ptf(l + T — 



then 

ll*2(t,*o,y 0 ) - *i(t,* 0 ,y 0 )ll < 
l|xi(t,* 0 .yo) -* 0 II + 



+ 



\\yi(t,x 0 ,y 0 ) -y 0 \\ 
And by mathematic induction we get : 



. . . (21) 



\ JL m + l 

\x m 



\\y m (t,x 0 ,y 0 ) -y m _i(t,x 0 ,y 0 )|| , 
and 



Hy2(t.*o-y 0 ) -yi(t,x 0l y 0 )\\ < 

\\xi(t,x 0 ,y 0 ) - x 0 \\ + 

+ 



t ATH 2 T\ i AT \ 

ATH 2 T\ ( AT \ 



(22) 

72-1 



RTF T\ ( BT \ 

FT -TT^) L i + 2pF ( 1 + TT^) L 

( RTF 2 T\ ( RT \ 

(Fr-^jL 2 + 2pF(l + TT ^)L 4 



\\yi(t,x 0 ,y 0 ) -y 0 || 

And by mathematic induction we get : 



lly m +i(t.*o.yo) -y m (t,x 0 ,y 0 )\\ < 
lk m (t^o.y 0 ) - ^ m -i(t.^y 0 )ll + 



. . . (23) 
ATH 2 T S 



( ATH T\ i AT \ 

[ HT -T^) K ^ 2 H 1 + T^h 

ATH 2 T\ i AT \ 

HT -1^) K2 + 2PH { 1 + 1^) K ' 



\\y m (t,x 0 ,y 0 ) -y m -i(t,x ol y 0 )\\ , ... (24) 
Rewrite inequalities (22) and (24) in vector form as : 
^ m +i(t^o.y 0 ) ^ A(t)V m (t,x 0 ,y 0 ) 
where 

/\\x m+1 (t, x 0 , y 0 ) - x m (t, x 0 , y 0 ) || N 

Vm+iit.x^yJ = 

\lly m +i(t^o.yo) -y m (t,x 0 ,y 0 )\\ / 

M(t) JV 3 (t) N 

A(t) = 

WO W 4 (t), 

and 

/||x m (t,x 0 ,y 0 ) -x m _!(t, x 0 ,y 0 )|P 

Kn(t.^o.y 0 ) = 

\lly m (t^o.y 0 ) -y m -i(t^o.y 0 ) 

It follows form the inequality (25) that : 
V m+1 < A 0 (t)V m 
where 

A 0 = max|A(t)| 

t£[0,T] 

By iterating the inequality (3.20) , we find that 
V m+1 <A™V 0 
where 



(25) 



. . (26) 
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// ATH T\ i AT \ 

( BTF 2 T\ ( RT \ 

Vr-rT^j M2 + 2pF ( 1 + iT^) M V 



which leads to the estimate 

m m 

^7i<^ArV 0 ...(27) 

1=1 1=1 

Since the matrix A 0 has Eigen -values like as (12), then the series (27) is uniformly convergent, i.e. 

m oo 

lim Y A 1 - 1 V 0 =Y A 1 - 1 V 0 = (E- . . . (28) 

m->cc / i / i 

i=l i=l 

where E is a unity matrix. 

The limiting relation (28) signifies a uniform convergent of the sequence x m (t,x 0 ,y 0 )andy m (t,x 0 ,y 0 ) in the 

domain (14). 

Let 

lim x m (t, x 0 , y 0 ) = x x (t,x 0 ,y 0 ) 



lim y m (t, x 0 , y 0 ) = y x (t, x 0 , y 0 ) 



(29) 



Since the sequences of functions x m (t, x 0 ,y D )andy m (t, x 0 ,y Q ) are periodic in t with period T , then the 
limiting of them are also periodic in t with period T , end thus x x (t, x 0 ,y 0 ) = x(t, x 0 ,y 0 ), y 00 (t, x 0 ,y 0 ) = 

y(t^ 0 .y 0 )- 

Also from (29) the following inequality 
f \\x x (t,x 0 ,y 0 ) -x m (t,x 0 ,y 0 )\\\ 

UAT^-^rt , ... (30) 

Jlya=(t,*o-y 0 ) -y m (t^o-yo)ll/ 

is hold for m > 1 , and hence 

*o=(t,*o.y 0 ) = *(t>*o.yo) . y»(t.^o.y 0 ) = y(t.x 0 .y 0 ) 

which are the solutions of the system (1). 



Uniqueness of Solution of (1) 

Let all assumptions and conditions of theorem 1 were given, then the two functions x(t,x 0 ,y 0 ) and 
y(t,x 0 ,y 0 )are uniqueness of solution of (1) in the domain (14). 
Proof : 
Let 

t 

x(t,x 0l y 0 ) = x 0 + j e Ht - s) [f(s,x(s,x 0 ,y 0 ),y(s,x 0 ,y 0 ))- 
o 

T 

r ^ 

- J 1 ~ glT e^ T - s) f(s,x(s,x 0 ,y 0 ),y(s,x ol y 0 ))ds]ds + 
o 

+ ^ e A(r - s) / i (x(t,x 0 ,y 0 ),y(t,x 0 ,y 0 )) - 

0<t ( <t 

p 

A \— ' 

~ i + e ^ 2_, e X(J ~ s) h (*(t, s 0 , y„), y (t, s 0 , y 0 )) , ... (31) 



i=l 

and 

t 



y(t,* 0 .y 0 )= yo + | e/?(t s) Lf( 5 ^( 5 ^o-yo).y(^s 0 ,y 0 )) - 

0 

r 

~ / 1 + eP T eHT ~ s) f( s >x( s > x °>y°)>y( s > x °>y°)) ds ] ds + 
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+ 



^ e p(J s} Gi(x(t,x 0 ,y 0 ),y(t,x 0 ,y 0 ))- 



0<tj<t 

- Y^wL e ~ s)G i{^ x o-yo)-yit.x 0 ,y 0 )) , ... (32) 
i=i 

Are another solutions for the system (1), then we shall prove that x(t, x Q , y 0 ) = x(t, x 0 ,y 0 ) , (t, x 0 ,y 0 ) = 
y{t, x ol y 0 ) , and to do this we need to prove the following inequality by induction, 
'||x(t,x 0 ,y 0 ) - x m (t,x 0 ,y 0 )\\\ /\\x(t,x 0 ,y 0 ) - x m - 1 (t,x 0 ,y 0 )\\\ 

)<K[ ■■■ (33) 

J\Kt,x 0 ,y 0 )-y m (t,x 0 ,y 0 )\\J \\\Kt,x 0 ,y 0 ) - y m -i(t,x 0 ,y 0 )\\J 
for m > 1 , where 

M i=r , m n ax n ll/(t.*-y)ll - M 2 = max ||flf(t,x,y)|| 



te[0,r] 



te[0,r] 



M 3 * = max n !!/;(*, y) II , M|= max NG^y) II 

l<i<p l<i<p 

For m = 0 in (15) and (16), we have 

ATH 2 T 



( ATH 2 T\ ( AT \ 

\\x(t,x 0 ,y 0 ) -x 0 \\< \HT - Y^W) M i* + 2 P H V 1 + TTi^J Ms * 

and 

/ RTF 2 T\ ( RT \ 

\m,x 0 ,y 0 ) - y 0 \\ < [FT - Y^pfj M 2 * + 2pF (1 + Y^ffr) M 4 

/ ATH 2 T\ i AT \ 

t ATH 2 T\ i AT \ 

\ HT -1^) K2 + 2 H 1 + 1^) K - 
■ (36) 



and for m = 1, we get also 
\\x(t,x 0 ,y 0 ) -xi(t,x 0 ,y 0 )|| < 
ll*(t,x 0 ,y 0 ) — II + 



(34) 
(35) 



lly(t,* 0 .yo) -y 0 ll 

and 

lly(t^ 0 .y 0 ) -yi(t.^y 0 )ll ^ 

\\x(t,x 0 ,y 0 ) - x 0 \\ + 



( PTF 2 T\ ( (3T \ 

/ RTF 2 T\ I RT \ 

{ FT -TT^) L i + 2pF ( 1 + TT^) L 



lly(t^ 0 .y 0 ) -y 0 ll 

Suppose that (33) is true for m = p -1 , i.e. 



. . . (37) 

|£(t<*o<yo) ~~ x p-z(.t> Xo'Vo) 



\x{t,x 0l y 0 ) — x p _ 1 (t, x 0 ,y 0 ) 

•<Ari 

K \\y{t,x 0l y 0 ) - y v _ 1 {t,x 0l y 0 )\\) \\\y(t,x ol y 0 ) - y v _ 2 {t,x ol y 0 ) 

then 



(38) 



\\x(t,x 0 ,y 0 ) - x p (t,x 0 ,y 0 )\\ < 

||*(t,*o.y 0 ) -^ P -i(t^o.y 0 )|| + 

+ 

||y(t,x 0 ,y 0 ) - y p _i(t,x 0 ,y 0 )|| 
and 



/ ATH 2 T\ ( AT \ I 

ATH 2 T\ i AT \ 1 



\\9(t, x 0 , y 0 ) - y p (t, x 0> y 0 ) || < (ft - L 1 + 2 P f(i + J?-^ L 

V RTF 2 T\ ( RT \ 

[ FT -YT^) L ^ 2pF ( 1 + YT^) L 



\x(t,x ol y 0 ) -x p _i(t,x 0 ,y 0 ) + 
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\\y(t,x 0 ,y 0 ) -y p -i(t,x 0 ,y 0 )\\ 
i.e. 

'\\x(t,x 0 ,y 0 ) - x p (t,x 0 ,y 0 )\\\ /\\x(t,x 0 ,y 0 ) - x p _i(t,x 0 ,y 0 )|| 

J|y(t.*o.y 0 ) -y P (t,x 0 ,y 0 )\\J \||y(t,x 0 ,y 0 ) -y p _i(t,x 0 ,y 0 )|| 

Then the inequality (33) is true for m = 0, 1, 2, . . . 
By iterating the inequality (33) gives : 

|£(t.* 0 .y 0 ) -* m (t.*o.yo)ll\ /P(t-* 0 >y 0 ) -* m -i(t.*o.y 0 ) 
Jly(t,*o,y 0 ) -y m (t.*o-yo)ll/ \lly(t,*o,yo) -y m -i(t^o.y 0 ) 

But from the condition (11) we obtain A™ -» 0 as m -» oo, hence, proceeding in the last inequality to the limit 
we obtain that x(t,x 0 ,y 0 ) — £(t,x 0 ,y 0 ) and y(t,x 0 ,y 0 ) = y(t,x 0 ,y 0 ) which proves that the two solutions 
x(t,x 0 ,y 0 )and y(t, x 0 ,y 0 ~) are unique, and this completes the proof of theorem 2. 

Existence of solution of (1) 

The problem of existence solution of the system (1) is uniquely connected with the existence of zero of 
the function A(0, x 0 , y 0 ) and A* (0, x 0 , y 0 ), which has the form: 

T 

A 



A f 

A(0,x o ,y o ) = 1 + e xr [ J e A(T - s) f(t,x trj (t,x 0 ,y 0 '),y m (t,x 0 ,y 0 ))dt + 
o 

p 

+ Y, e ^ T ~ s)I ^ x ^ t i- x o-yo^-y^ t i- x o-yo))] ■ ■ ■ ( 39 ) 

i=i 

T 

A*(0,s„,y o ) = 1 / e/? r [ / e^^/^^^^yj.y^^^yj)^ + 



o 

p 



+ ^e' ? ( r - 5 )Gi(x 00 (ti,x 0 ,y 0 ),y 00 (ti,x 0 ,y 0 ))] . . . (40) 

i=i 

Where x x (t, x o ,y o )andy O0 (t, x D ,y 0 ) the sequence's limiting (12) and (13) successively, and this function is 
approximately determined from the sequence of functions: 

T 

A r 

A m (o,* 0 .yo) = 1 - e „ [ J e 1(7, - s) /(t^ m (t^o.yo).y m (t^o'yo))rft + 

0 

p 

+ ^e^ r - s )/ i (x m (t i ,x 0 ,y 0 ),y m (t i ,x 0 ,y 0 ))] . . . (41) 

;=i 

r 

Am(Q.^o.yo) = 1 f e pr [ / e/?(r " s) /(t^ m (t^o-yo).y m (t^o-y 0 ))rft + 



0 

p 



+ ^e^ r -^G i (x m (t i ,x 0 ,y 0 ),y m (t i ,x 0 ,y 0 ))] . . . (42) 

i=i 

where m — 0,1,2,... 
Theorem 3 : 

Let all assumptions and conditions of theorem 1 were given, then the following inequality : 
'|| A(0,x o ,y o )-A m (0,x o ,y o ) 

< Q A™(£ - A 0 ) -1 l£ ...(43) 

||A*(t,x o ,y o )-A^(0,x o ,y o 
where 

/ AHT ApH AHT ApH \ 

7^7 Ki + — Ko — Kn + TZtKa 

l + e AT l + e AT l + e AT 1 + e AT 

PFT PpF fiFT PpF 

Is holds for m > 0 , t £ [0, T] , x 0 £ % , y 0 £ % 
proof : 
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According to (39) and (41) , we have. 



A C 

l|A(0,x o ,y o )-A m (0,x o ,y o )|| < j IK^II \\f(t,x x (t,x 0 ,y 0 ), yoo (t,x 0 ,y 0 )) - 



-fit, x m (t, x 0 , y 0 ), y m (t, x ol y 0 ))\\dt + 



- h (x m {t u x 0 , y 0 ), y m {t u x 0 , y 0 )) || < 
AHT 



A \— ' 

+ YVe"lJ eKT ~ S) " " 7i ^°° (ti ' X °' yo) ' y °° (ti ' X °' ~ 
i=i 



+ 



l + e lT 
ApH 



[K^x^t.x^y^ - x m (t,x 0 ,y 0 )|| + K 2 \\y x {t,x 0l y 0 ) - y^t.x^y^W + 



So that 



l + e 1T 



■ [^3 Ikoo (t- * 0 . Jo) - *m 0> * 0 . y 0 ) II + K 4 \\y x (t, x 0 , y G ) - y m (t, x 0 , y 0 ) ||] 

II A(0,x o ,y o ) - A m (0,x o ,y o )|| < {^^W K i + fX^f ^) H^ft^X)) - *m (t. *<>. X>) II 

/ AWT Aptf \ 

+ u + e AT Kl + 1 + e 17, 1(4 ) Wy^'Xo-yo) -ym(t,x 0 ,y 0 )\\ 



(44) 



By the same method and by (40) and (42) , we have 

||A*(t,x 0 ,y 0 ) - A^(0,x o ,y o || < - ^ L x + - - ^ L 3 J ^xJs,x 0 ,y^ - x m {t,x 0 ,y 0 )\\ 

( PFT fi-pF \ 

+ u + eP T Ll + 1 4- eP T 14 ) \\y™(t- x o-yo) - ym(t.x 0 ,yo)\\ 

... (45) 

And so on, rewrite the inequalities (44) and (45) in vector form as : 



<Q 



\\x x (t,x 0 ,y 0 ) - x m {t,x 0l y 0 )\\ 
,\\ya,(t,x 0 ,y 0 ) - y m (t,x o ,y 0 )\\ 



<QA™(E-A o y% 



'II A(0,x o ,y o )-A m (0,x o ,y o )|| 

l|A*(t.^o.yo)-^(0,x o ,y o || 
And by (30) , we get 
'|| A(0,x o ,y o )-A m (0,x o ,y o )|| N 

||A*(t,x o ,y o )-A^(0,x o ,yJ 
Theorem 4 : 

Let thefunction f(t,x,y) and g(t,x,y) in the system (1) are defined on the intervals [a, b]and [c, d] 
respectively, and periodic in t with period T. 

Let that the sequence of functions (41) satisfying the next inequalities : 
min A m (0,x o ,y o ) < -S m \ 
a + M 5 <x 0 <b-M s \ 

max A m (0,x o ,y o ) > S m 
a + M 5 < x 0 < b - M 5 
Let that the sequence of functions (42) satisfying the next inequalities : 
min A^(0,x o ,y o ) < -e m ^ 
c + M 6 <y 0 <d-M 6 



(46) 



max A^(0,x o ,y o ) > e m 
c + M 6 <y 0 <d-M 6 

for m > 0 where : 

AHT ApH \ 1 AHT 



(47) 



KAHT ApH \ I AHT ApH \1 „ 

_ W PFT t PpF 



RFT RpF 

■L, + . „„ L, 



a + eH* " 2 " 1 + el> T UA J 

Then the system (1) has periodic solution of period T, x — x(t, x 0 ,y D )and y = y(t, x 0 ,y D )for which a + M 5 < 

x 0 <b-M 5 , c + M 6 <y 0 <d-M 6 . 

Proof: 

Let X\ , x 2 be any two points in the interval [a + M Sl b — M 5 ] such that : 
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A m (0,*i,yi) = min A m (0,x o ,y o )^ 
a + M 5 < x 0 < b - M 5 

(48) 

(0,x 2 ,y 2 ) — max A m 
a + M 5 < x 0 < & - M 5 
Let yi , y 2 be any two points in the interval [c + M 6 , d — M 6 ] such that : 
A* m (0,x 1 ,y 1 ) = min A^ (0, x 0 , y 0 ) ^ 
c + M 6 <y 0 <d-M 6 

A^(0,x 2 ,y 2 ) = max A^(0,x o ,y o ) 
c + M 6 <y 0 <d-M 6 j 

By (43) and (46), we get : 

A(0,*i,yi) = A m (0,x 1 ,y 1 ') + [A(0,x 1 ,y 1 ) - A m (0,x 1 ,y 1 )] < 0 

A(0,x 2 ,y 2 ) = A m (0,x 2 ,y 2 ) + [A(0,x 2 ,y 2 ) - A m (0,x 2 ,y 2 )] > 0 
From (43) and (47), we get : 

A*(0,x 1>yi ) = A* m {Q,x 1>y {) + [A\0, Xl , yi ) - A* m {Q,x 1>yi )] < 0 



■ ■ ■ (49) 



(50) 



. . . (51) 



A*(0,x 2 ,y 2 ) = A;(0,x 2 ,y 2 ) + [A*(0,x 2 ,y 2 ) - A^(0,x 2 ,y 2 )] > 0 
It follows from the functions (39) , (40) and the relations (50) , (51) in virtue of the continuity of the A- 
constant, that there exists x x — x Q , x x £ [%,x 2 ] and y x — y 0 , y x £ [yi,y 2 ] such that A(0,x oo ,y 0O ) = 0 , 
A*(0,x a3 ,y OD ) = 0. And this proved that the system (1) has a periodic solution x = x(t,x 0 ,y 0 )for x 0 £ [a + 
M 5 ,fo-M 5 ]andy = y(t,x 0 ,y 0 )fory 0 £ [c + M 6 ,d-M 6 ]. 
Remark 1 161 : 

When R n = R , i.e. when x is a scalar theorem 4can be strengthens by giving up the requirement that 
the singular point shout be isolated, thus we have 
Theorem 5 : 

Let the function A(0,x o ,y o ) defined asA : D X f -» R 

T 

X 



A(0,x o ,y o ) = 1 + e xr i J e l( - T - s) f(t,x m {t,x 0 ,y 0 ),y x {t,x 0 ,y 0 ))dt + 
o 

p 

+ ^e^ T ~ s: >I i (x 00 (t i ,x 0 ,y 0 ),y 00 (t i ,x 0 ,y 0 ))] 

i=l 

Let the function A*(0, x ol y 0 ) defined asA* : Dp f -» R ' 

T 

A*(0,x o ,y o ) = 1 + e p T [ j e^ T ~ s ' ) f(t,x aD (t,x 0 ,y 0 '),y aD (t,x 0 ,y 0 ))dt + 
o 

p 

+ X eP iT ~ S) Gi ^' X °' y °*' y ™ ^ i,Xo ' 



(52) 



(53) 



Where the functions x a) (t,x 0 ,y 0 )andy 0= (t, x 0 ,y 0 )are the limit of asequence of periodic functions (12) , 

(13) respectively, then the following inequalities are holds : 

/ AtH 2 T\ ( At \ 

l|A(0,x o ,y o )|| <(//t- T - r ^ F jM 1 +2p//^ + T - r ^ F jM 3 . . . (54) 

/ ptF 2 T\ ( fit \ 

HA*(0,x o ,y o )|| < [Ft - M 2 + 2pF (l + ^-^J M 4 



(55) 



UAOU^yi) - A(0,x 2 ,y 2 )|| < 1^(1 - W±- W 2 W 3 (1 - I^)" 1 ) _1 [\\x$ - x 2 \\ + 

+w 2 (i-iy 4 )- 1 ||y 0 1 - y < 2 ||] + 
+w 2 {i-w ± -w 2 w 2 {i - w.y 1 r 1 [llyo 1 -yl II + 
+WzO- - w.y 1 ^- xl ||]_ . . . (56) 

HA'CO.^.yi) - A*(0,x 2 ,y 2 )|| <W,(1 -W, - W 2 W 3 (1 - M4) _1 y 1 [\\xl - x 2 0 \\ + 

+W 2 (l-W i y 1 \\yi- yZ\\] + 
+VU1 -W4- W 2 W 3 (1 - Wi)" 1 y 1 [Hyo 1 - y 2 || + 
+W 3 (l-W 1 )- 1 ||x 0 1 -x 2 ||] ...(57) 

Where 
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W 2 



\( ATK X \ ( AK 3 p \ 



AK 3 p 



AK 4 p 



W 3 = FT 



W 4 = FT 



( l i-it^) + ( L3+ itS^) 



PTL 2 



2 l + ePT 
for all x 0 , ,xl £ % and y 0 ,yl ,y 2 £ %. 

Proof: 



From the properties of the function^ (t, x 0 ,y 0 )andy 0G (t, x Q , y 0 ) established by theorem 1, it follows that the 
functionsA= A(0, x 0 , y 0 ) and A*= A* (0, x 0 , y 0 )are continuous and bounded with the positive constants 
/ AtH 2 T\ ( At \ 

\ Ht ~ TTe^j Ml + 2pH V + TT^J M3forx ° £ Dlf • and 

[Ft - 1 + efST ) M 2 + 2pF (l + 1 + efST ) M 4 fory 0 £ % , respectively. 
By using (52) , we have 

T 

A f 

UAOU^yi) - A(0,x 2 ,y 2 )|| < J \\e^ T -^\\\\f{t,x x {t,xl,yl),y x (t,xl,yl)) ~ 

o 

-f{t, x x (t, x 2 , y 2 ), y x (t, x 2 , y 2 )) \\dt + 

A P 

+ 1 + elr ^||e A(r - s) ||||/ i (x 0O fe,x o 1 ,y o 1 ),y O3 (t i ,x o 2 ,y o 2 )) - 
i=i 

- k (x^, xl,yl), y ao {t i ,xl, y 2 )) || < 
AHT 

< [Ki\\x a> (t,4,yZ) - Xoo(t.* 0 2 .yo)ll + K 2 \\y x (t,xl,yl) - y x (t,x 2 ,y 2 )\\] + 

1 + e Al 

+ -^^VK 3 \\x m {t,xl,yl) -x x (t,x 2 0 ,y 2 )\\ + K A \\y x (t,xl,yl) - y x {t,x 2 ol y 2 )\\\ 



1 + e 



So that 



AH 



UACO^yi) - A(0,x 2 ,y 2 )|| < ^—^[(77^ + pK 3 )\\x m (t,xly^ - x x {t,x 2 0 ,y 2 )\\ + 

HT^+pK^Wy^t.xly^-y^t.xly 2 )]]] . . . (58) 
And we will find by the same method and by (53) , we have : 

RF 

IIA'CO.Xo 1 , y,, 1 ) - A*(0,x 2 ,y 2 )|| < -^—^ [QTL, + pL^Wx^t.xly^ - x x {t,x 2 0 ,y 2 )\\ + 

+ (rL 2 +pL 4 )||y OD (t,x o 1 ,y o 1 )-y O0 (t,x 2 ,y 2 )||] . . . (59) 
Where x ca (t, xl, yl ) , x oc (t,x 2 ,y 2 )andy 30 (t,x 0 1 ,y 0 1 ) , y OD (t,x 2 ,y 2 ) are the solutions of the following integral 
equations : 



x(t,x 0 fc ,y 0 fc ) = *o + j" e Ht s) [f(s,x(s,xtyo),y(s,x*,y*)) - 



o 

T 



__ e HT ^f(s,x(s,x*,y*),y(s,x*,y*))ds]ds + 



o 

+ £ eW-^I^xiU.xlyUyiU.xly*))- 



o<t,<t 



1 + e 



e^^k (x^xly*), y(t ; , x k a , y 0 k )) 



(60) 



And 



y(t,x k ,y k )= y 0 fc + | e' 3(t - s) [/(5,x(s,^,y 0 ' c ),y(5,x 0 fc ,y 0 ' c )) - 



i 

-/ 



l + eP T 



e P{T s) f(s iX ( SiX k ,y o fc ),y( 5 ,^ ,y*))ds]ds + 
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+ 

0<tj<t 

P 



l + e? 7 

Where k = 1,2 
From (60) , we have 



^e^^G^x^xly^.y^xly*)) , ...(61) 



\\x a> (f,xl,yl)-x x (f,xl,yl)\\= xl + J e A(t s) [/(5,x O0 (5,x o 1 ,y o 1 ),y 3O (5,x o 1 ,y o 1 )) - 



/ 



T 

1(r - s) /(s,x 0O (5,x o 1 ,y o 1 ), y^O,**, y D 1 ))ds]d5 + 



l + e Ar 

o 



+ 2^ e^-^liixMi.xl.yD.ySU.xl.yl))- 
p 

^ e ;i(r - s) /i(^(ti,^,yo 1 ).y»(c<.^.yo 1 )) - 



0<tj<t 

p 

1 + e^ 

i=l 

t 



*o-J e 1(t s) |/(s,x«,(s.*o.yo).y»(s»*o.yo 2 )) - 

0 

T 

f A 

J i + e Ar eA(r ~ 5) /(^^oo(s,^,y 0 2 ),y 00 (s,^,y 0 2 ))ds]ds - 
o 



^ e Ap, - s) / < (ac G0 (t( > ^ f y o 2 ) f y G0 (t( f x2 f y2)) + 



0<tj<t 

p 

A 



+ i + e Ar X eA(r ~ 5)/ '(*°°fe'*o.yo 2 ).yoo(t,x 2 ,y 0 2 )) < 

i=l 



< ll^ 1 - x 2 || + ///f^ llacofox^yo 1 ) - ^(t.^.y 2 )!! + 

At/zr/fi 

+ ,~ , Tr \\y«>(.t,xlyZ) -yoo(t,x 2 ,y 2 )|| + tf/f 3 t ||* x (t,^ > y 0 1 ) - x Ta {t,xl,y^)\\ + 



+ //tf 2 t \\y x (t,xl,yl) - y„(t,x 2 ,y 2 )|| + lr Ux^t^y, 1 ) - y 0 2 )|| + 



AtHTK 2 
1 + e 1 ^ 



+ ///f 4 t Hy^t^y, 1 ) -y„o(t,x 2 ,y 2 )|| + ^^-p\\x x (t,xl,y^) - x x (t,xl,y*)\\ + 

1 + e M 



+ ^f r fr PllyooCt^o-yo 1 ) - yoc(t,x 0 2 ,y 0 2 )ll 

1 + e Al 

< - x 2 || + tft ^ + 1 + elr +K 3 + 1 + eAT ) \\x„(t,xl.yZ) ~ ^(t,xly^)\\ + 

+ //t (k 2 + + /f 4 + y^w) HyooCt.^.yo 1 ) - ySt,x 2 0 ,yt)\\ . . . (62) 

By the same method and by (61), we find : 
\\y x (t,xl,yl) - y x (t,xl,yt)\\ < \\yl -y 2 || + 

+ Ft (l x + ^ L ^ T + L 3 + ^f pT ) WxSt.xl.yl) - x x {t,xl,y*)\\ + 

+ Ft (L 2 +^^ + L4+f^7jlly»(t^o 1 ,yo)-y»(t,x 2 ,y 2 )|| . . . (63) 
From (62), we have: 

, , / / AT/C, Ar/fopxA -1 
Ikoo(t^o.yo) - ^(t,x 2 ,y 2 )|| < 1 1 - //t (#! + ^-^r + /f 3 + Y^rjw) J 

(ll^o - *o 2 II + m (/c 2 + + ^ + lly»(t.*o.yo) - y»(t,* 2 ,y 2 )||) 

1 ... (64) 

By (63), we have: 
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Ily- it, xl yl) - y^t.xlyl) || < (l -Ft (L 2 + I ^-^ + L 4 +^^)j 

(llyo 1 - yl II + f t (l x + + l 3 + ^t^f) Hy- d. *h yi) - y*> * 2 . y 0 2 ) ll) 

... (65) 

By substitute (65) in (64) , we obtain :W 3 (1 - W 4 ) _1 

ll^ooCt^o.yo^-^ooCt^o.yo)!! ^ - ^\ + W 1 \\x 0 ,{t,xl,yl) - x^it.xl.yDW + 
+ W 2 (l - ^ 4 )- 1 (lly 0 1 -y 0 2 ll + WzWx^t.xl.yl) - x a {t,xl,yS)\\) 

So that 

\\x x {t,xl,yl) - x x {t,xl,yl)\\ < (1_- W x -WffiQ. - W^YWUl - x 2 || + 

+w 2 (i-w 4 )- 1 ||y 0 1 -y 0 2 ll) ...(66) 

And also, we have 

HyooCt^o'.yo^-ycM.y 2 )!! < \\yl -yjjl + w~ 4 \\y a ,(.t,xlyt) -y a (t,xly^\\ + 

+ W 3 (1- Wj^Qxl - x 2 0 \\ + W 2 \\y O0 (t,x 1 0 ,y^ - y M (t,x 2 ,y 2 )||) 

So that 

\\yUt,xlyZ)-y m (t,xlyZ)\\ < (l- -WMa - TPD" 1 )- 1 ^ 1 - y 2 ll + 

+W 3 (l-Wi)- 1 ||x 0 1 -x 0 2 ||) ...(67) 
Now by substitute (66) and (67) in (58) and (59), we obtain on (56) and (57) respectively. 
Remark 2 Til: 

The theorem 5 to ensure solution's to the system (1), in view of to happen small change in the point x 0 
, to requite small change on the function's behaviorA= A(t,x 0 ,y 0 ). 
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